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Abstract 

We investigate magnetic properties and effects of pairing fluctuations in the BCS (Bardeen- 
Cooper-Schrieffer)-BEC (Bose-Einstein condensation) crossover regime of an ultracold Fermi gas. 
Recently, Liu and Hu, and Parish, pointed out that the strong-coupling theory developed by 
Nozieres and Schmitt-Rink (NSR), which has been extensively used to successfully clarify various 



physical properties of cold Fermi gases, unphysically gives negative spin susceptibility in the BCS- 
BEC crossover region. The same problem is found to also exist in the ordinary non-self-consistent 
T-matrix approximation. In this paper, we clarify that this serious problem comes from incom- 
plete treatment in term of pseudogap phenomena originating from strong pairing fluctuations, as 
well as effects of spin fluctuations on the spin susceptibility. Including these two key issues, we 
construct an extended T-matrix theory which can overcome this problem. The resulting positive 
spin susceptibility agrees well with the recent experiment on a 6 Li Fermi gas done by Sanner and 
co-workers. We also apply our theory to a polarized Fermi gas to examine the superfluid phase 
transition temperature T c , as a function of the polarization rate. Since the spin susceptibility is 
an important physical quantity, especially in singlet Fermi superfluids, our results would be useful 
in considering how singlet pairs appear above and below T c in the BCS-BEC crossover regime of 
cold Fermi gases. 

PACS numbers: 03.75.Ss, 03.75.-b, 03.70.+k 
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I. INTRODUCTION 



The uniform spin susceptibility x is a fundamental quantity in considering magnetic 
properties of an electron system. In a free electron gas, \ gives useful information about 
the single-particle density of states at the Fermi level In s-wave superconductivity, x is 
suppressed below the superconducting phase transition temperature T c to vanish at T = 

, because the spin degrees of freedom become inactive by the formation of singlet Cooper 
pairs. The suppression of the spin susceptibility has been also observed in the underdoped 
regime of high-T c cuprates, which is referred to as the spin gap phenomenon in the literature 

• Although the origin of the spin gap is still in debate, the importance of preformed pairs 



has been pointed out 
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-8| Fermi gases, the high tunability of 
Indeed, using a tunable pairing 



Since the realization of superfluid 40 K 5] and 6 Li 
this quantum system has attracted much attention 

interaction associated with a Feshbach resonance [12|, one can study superfluid properties 
from the weak-coupling BCS regime to the strong-coupling BEC limit in a unified manner 

13)- 17]. In the so-called crossover region, a deviation of single-particle 



(BCS-BEC crossover) 



excitation spectrum from the free particle dispersion has been observed in the normal state, 
by using the photoemission-type experiment developed by JILA group [18|, |l9|. As an 



explanation for this anomaly, the possibility of the pseudogap phenomenon associated with 



strong pairing fluctuations has been proposed |20l-|27j. Since the cold Fermi gas system is 
much simpler than high-T c cuprates, the former system would be useful for the assessment 
of the preformed-pair scenario discussed in the latter. 

Besides the tunable interaction, the high tunability of population imbalance is another 
advantage of cold Fermi gases {28, 



29]. When we describe two atomic hyperfine states in 



a Fermi gas by pseudospin a =t, I, a polarized Fermi gas is closely related to an electron 
system under an external magnetic field. In the limit of low population imbalance, one 
may evaluate the spin susceptibility. Using this quantity, one can examine whether the 
preformed singlet pairs really appear in the BCS-BEC crossover regime of a cold Fermi gas. 
In the case of a finite population imbalance, the mismatch of the Fermi surfaces between 
the t _s P m component and |-spin component is expected to cause the instability of the 



s-wave superfluid state 



30|, where various exotic states have been proposed, such as the 



Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state 



31 



32| and the Sarma phase [33M35 ] . 
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In this paper, we investigate (pseudo)magnetic properties of a normal state Fermi gas 



in the BCS-BEC crossover region. In the un po 
developed by Nozieres and Schmitt-Rink (NSR) [15 
clarify various physical properties of this system [17|, 



arized case, the strong-coupling theory 



ras been extensively used to successfully 
2o|, bli bd, U, bd-kll • However, when 
we apply this theory to a polarized Fermi gas, it is known that negative spin susceptibility 



is obtained in the crossover region 



42 



431 ] (which is thermodynamically forbidden |44|). 
Because of this serious problem, so far, the phase diagram of a polarized Fermi gas has 
mainly been examined within the mean-field level 35[. However, as in the unpolarized 
case, strong-coupling effects would be also important in the BCS-BEC crossover regime of a 
polarized Fermi gas. Indeed, it has been pointed out that the FFLO state (which has been 



predicted in a polarized Fermi gas within the mean-field analysis 
pairing fluctuations 
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351 ] ) is unstable against 



46] . Thus, to discuss the BCS-BEC crossover physics of a polarized 



Fermi gas, we need a reliable and tractable strong-coupling theory which can overcome the 
above mentioned problem. 

In this paper, we show that the "negative susceptibility problem" also exists in the 
ordinary (non-self-consistent) T-matrix approximation, which has been also extensively used 
in the unpolarized case. Clarifying the origin of this serious problem, we present a minimal 
extension of the T-matrix theory to correctly give the required positive spin susceptibility 
in the whole BCS-BEC crossover region. The calculated spin susceptibility in this extended 
T-matrix theory is shown to agree well with the recent experiment on a 6 Li Fermi gas j^ . 
We also apply this theory to the system with finite population imbalance, and examine the 
critical population imbalance at which the superfluid phase transition disappears. 

This paper is organized as follows. In Sec. II, we explain our formulation. We also 
compare our theory with the NSR theory, as well as the ordinary T-matrix approximation. 
In Sec. Ill, we calculate the spin susceptibility to show that our strong-coupling theory does 
not meet the negative susceptibility problem in the whole BCS-BEC crossover region. We 
also compare our results with the recent experiment on a 6 Li Fermi gas. In Sec. IV, we treat 
a polarized Fermi gas. Throughout this paper, we set H = &b = 1, and the system volume 
V is taken to be unity. 
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II. MODEL POLARIZED FERMI GAS AND STRONG-COUPLING THEORIES 



We consider a two- component Fermi gas with population imbalance, described by the 
BCS Hamiltonian, 

H = £p,v C h° C P,°- ~ U ^2 C l+Q/2,t C ~p+q/2,i C -p'+Q/2,i C p'+q/2^- (!) 

p,a P,p',q 

Here, c pa is the creation operator of a Fermi atom with momentum p and pseudospin 

2 

a =t, 4, describing two atomic hyperfine states. £ P)CT = e p — fi a = #- — fj, a is the kinetic 
energy of the cr-spin component, measured from the Fermi chemical potential \x a (where m is 
an atomic mass). The pairing interaction —U (< 0) is assumed to be tunable by a Feshbach 
resonance. As usual, we measure the interaction strength in terms of the s-wave scattering 
length a s , given by 

47TO, -U 

where Co> c is a high-energy cutoff. In this scale, the weak-coupling BCS regime and the strong- 
coupling BEC regime are characterized by {kpCbs)^ 1 ^ — 1 and (/ci?a s ) _1 > 1, respectively. 
(Here, kp = p^iV] 1 / 3 is the Fermi momentum, where N is the total number of Fermi 
atoms.) The region —1 < {k F a s )~ l < 1 is called the crossover region. In this paper, we 
consider a uniform Fermi gas, for simplicity. 

When we write the chemical potential \i a as fi a = // + ah [where // = (/i^- + /xj.)/2 is 
the averaged chemical potential], Eq. (Q]) may be viewed as a model Hamiltonian for an 
interacting electron system under an external magnetic field h. The spin susceptibility x is 
then given by 

X = lim T - l . 3 

ft->o h 

Here, N a is the number of Fermi atoms with cr-spin, which is calculated from the single- 
particle thermal Green's function G Pta (iu n ) as 

N a = Tj2 G P A^n), (4) 

p,iu} n 

where u n is the fermion Matsubara frequency. In this formalism, strong-coupling effects on 
X is described by the self-energy S p cr (iu; n ) in G Pta (icu n ), 

[G pya {lU n )\ - £p 1<T (^n) 



(a) 



^*^% G q _ p _ {iv n -im n ) 



(b) 



► s ► K-E-j* 

Af>„) = \-U+ > > + 



FIG. 1: (a) Self-energy £ P)(r (zu; n ) used in this paper, (b) Particle-particle vertex function T(q,iu n ). 
The solid line and the solid double line represent the free Green's function G° and the full Green's 
function G in Eq. respectively. The wavy line describes the attractive interaction —U. 



Here, G pa {iuj n ) = [uo n — £ p ,o-] _1 is the Green's function for a free Fermi gas. 

As mentioned in the introduction, the NSR theory breaks down for a polarized Fermi gas 
in the sense that it incorrectly gives the negative spin susceptibility (x < 0) in the BCS- 



BEC crossover region |42|, |43| . This implies that one needs to carefully treat the self-energy 
correction Tj p ^(iuj n ) in considering magnetic properties of a polarized Fermi gas. In this 
paper, we take the strong-coupling corrections diagrammatically described by FigJT] (We 
will explain the reason for this choice in Sec. III.), which gives 

where v n is the boson Matsubara frequency. T(q, iv n ) is the particle-particle vertex function 
in the ladder approximation (See Figfjjb).), 

T(q,iu n )= , (7) 

1 - UU{q,iv n ) 



where 



U(q, iv n ) =T^ Gp+q/2,tOn + iuJ n)G°_ p+q/2il (-iu n ) 

p,iu)„ 



1 - /(£p+g/2,t) - f{tL-p+q/2,l) 



(8) 



W n — £ p +qr/2,t — i-p+q/2,\ 



is the lowest order pair propagator. In Eq. (151) , f(x) is the Fermi distribution function. 

The ordinary (non-self-consistent) T-matrix approximation (TMA) also uses the self- 
energy in FigJU except that the full Green's function G in FigJT^a) is replaced by the 
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noninteracting one G°, as 



(9) 



q,iv n 



In this sense, our strong-coupling theory may be regarded as an extended T-matrix approx- 



imation (ETMA) [48(| . We briefly note that, although the NSR theory also uses £p iCr i 



lUJr, 



the Green's function in Eq. (j5]) is expanded to 0(S°) as 



(10) 



As usual, the superfluid phase transition temperature T c is determined from the Thouless 
criterion, 

r- 1 (q,^ n = o)| T=:rc = o. (11) 

While the uniform superfluid state corresponds to q = 0, the FFLO state is realized when 
the highest T c is obtained at q ^ 0. However, since the latter is known to be unstable against 



pairing fluctuations even for a weak interaction in the absence of a optical lattice [45|, |46| , 
we set q = in Eq. ffTTj) from the beginning. In this case, the (regularized) T c -equation is 
given by 



m 

47T<2, 



tanh 



£p,t 
2T 



tanh ( % 
2T 



2e T 



(12) 



where ^ p = e p — is the kinetic energy, measured from the averaged chemical potential 
fi = (/if + /i|)/2. For a given total number of Fermi atoms N = + N±, we solve Eq. (|T2|) . 
together with the number equation (j3J), to determine T c , /x, and /i, self-consistently. In the 
unpolarized case, the three strong-coupling theories (ETMA, TMA, and NSR) qualitatively 
give the same BCS-BEC crossover behavior of T c , as shown in Figj2j In the next section, 
however, we show that they give very different results for the spin susceptibility. 



III. SPIN SUSCEPTIBILITY IN THE BCS-BEC CROSSOVER REGION 

Figure [3] shows the spin susceptibility x a t T c in the BCS-BEC crossover. As mentioned 
previously, the NSR theory gives the negative spin susceptibility (xnsr < 0), when the 
interaction becomes strong to some extent. The situation becomes better in the ordinary 
T-matrix theory (xtma)- However, as shown in the inset of Fig 01 Xtma slightly becomes 
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FIG. 2: (color online) Calculated T c of an unpolarized Fermi gas in the extended T-matrix approxi- 
mation (ETMA). For comparison, we also show the results in the ordinary T-matrix approximation 
(TMA), as well as the NSR theory (NSR). The interaction strength is measured in terms of the 
inverse scattering length a s , normalized by the Fermi momentum fcp. ep is the Fermi energy. 



negative in the crossover region. In contrast, our extended T-matrix approximation (xetma) 
gives the required positive spin susceptibility in the whole BCS-BEC crossover, xetma 
decreases with increasing the interaction strength, which reflects the increase of preformed 
Cooper pairs at T c . Since all the Fermi atoms form tightly bound singlet molecules in the 
BEC limit, Xetma vanishes in this limit. 

To understand the reason why the present ETMA can overcome the negative susceptibility 
problem, it is helpful to note that strong-coupling effects on \ can be divided into the self- 
energy part and the vertex part, as diagrammatically shown in FigH^a). Between the two, 
the former comes from the self-energy correction H p ^(iu n ) in the single-particle Green's 
function in Eq. (J5J), so that this part physically describes how strong-coupling effects on 
single-particle excitations affect the spin susceptibility \- 111 this regard, we recall that 
strong-pairing fluctuations cause the pseudogap phenomenon in the crossover region 20- 



27| . where a gap-like structure appears in the normal state density of states p(u) around 
the Fermi level u = 0. Since % is deeply related to p(0) 1^9 1 , the pseudogap leads to the 
suppression of x m the crossover region. 

However, as pointed out in Ref. |20 |. the NSR theory overestimates the pseudogap to 
incorrectly give the negative density of states around u = 0. This is because of the fact 
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FIG. 3: (color online) Calculated spin susceptibility x & t T c in the BCS-BEC crossover, xetma: 
extended T-matrix approximation, xtma: ordinary T-matrix approximation, xnsr: NSR theory. 
The asymptotic form of xnsr m Eq. (fl3j) in the NSR theory is also shown. The inset shows xtma 
magnified in the crossover region where it becomes negative, xo is the spin susceptibility of a free 
Fermi gas at T = 0. 



that the NSR theory only retains the self-energy correction to 0(S°). Thus, the NSR spin 
susceptibility Xnsr also becomes negative in the crossover region where the pseudogap be- 
comes remarkable in p(oS). Using the NSR Green's function in Eq. ([TO]) , one finds that Xnsr 
is diagrammatically gives by Figj4^b). In this panel, the second term (= Xnsr.) describes 



the pseudogap correction to % 



50] . which becomes dominant over the third term (which 



describes a vertex correction) in the BEC regime. In the BEC limit, one finds 



Xnsr — 



- , (b2) 

Xo + Xnsr 



167ra= 



Xo 



rn 



2mT BEc y / 3 
2^ J ^ \2 



d 2 N° 



dh 2 



(13) 



h=Q 



(We summarize the derivation in Appendix |A]) Here, 

Xo = — ~ sech 2 ( -2- | 
2T ^ \2T J 



(14) 



is the spin susceptibility of a non-interacting Fermi gas. In Eq. (TT3~j) . N® = J2 p f (£ p .t) is the 
number of j-spin atoms in a free Fermi gas. T C BEC = 0.218£f is T c in the BEC limit 15Nl7|. 
Since the non-interacting part Xo i n Eq- f|T3|) is remarkably suppressed in the BEC regime 



(a) 



(b) 

^nsr" 

(c) 
^TMA" 



OC> + <H> 



(bi) 



(b2) 



(b3) 
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FIG. 4: (a) Feynman diagram describing spin susceptibility x- The solid double line is the full 
Green's function involving the self-energy correction. A is a three-point vertex part, (b) xnsr- We 
only retain the terms to O(S ). (c) xtma- (d) Xetma- In panels (b)-(d), Y is the particle-particle 
scattering matrix in the ladder approximation in FigQTh). 



due to the negative chemical potential (/z < 0) |l4Hl7l|. the correction term Xnsr l ea ds to 
the negative spin susceptibility, as shown in Fig. [3j 



The pseudogap effect on p(u ~ 0) is correctly treated in TMA |20M27l|. However, this 
approximation still has a problem in the vertex part A, so that xtma becomes negative in the 
crossover region. To see the origin of this, we diagrammatically compare Xtma (where the 
self-energy E° in Eq. fl9]) is used) with Xetma (where the self-energy E in Eq. fl6]) is used) in 
FigJU While Xetma involves the random phase approximation (RPA)-like series of the Maki- 
Thompson (MT) diagrams 



50 



5l|, TMA only retains this series to the first order. When 
we approximate the particle-particle scattering matrix T to the bare interaction —U, and 
ignore all the other interaction effects, Xetma in FigH^d) reduces to the RPA susceptibility, 



Xetma 



Xo 



(15) 



1 + Uxo 

That is, the vertex part Aetma = + Uxo], as well as Xetma, are always positive. In 



contrast, because of 

Xtma ^ Xo[l - Ufa], 
Xtma becomes negative, when the vertex part A TM a = 1 



(16) 



Uxo becomes negative 52]. 
Since the present ETMA correctly treats both the self-energy part and the vertex part, the 
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FIG. 5: (color online) Calculated spin susceptibility \ hi the normal state (solid line). The exper- 
imental data [471 ] are shown as the solid circles. To reproduce the experimental situation [47J, the 
temperature is fixed at the value of T c for (fopa^) -1 = —0.35 (vertical line). While the left side of 
the vertical line is the normal phase, the right side is the superfluid phase. 
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FIG. 6: Calculated spin susceptibility x m the normal state above T c . We take (kFa s ) 1 = 0. Near 
T c , the decrease of % with decreasing the temperature is due to the pseudogap effect. 



required positive spin susceptibility is obtained over the entire BCS-BEC crossover region, 
as shown in Figj3j 

In Fig. |5l we compare the calculated spin susceptibility Xetma with the recent experiment 
on a 6 Li Fermi gas [43]. In this experiment, the temperature is fixed at the value of T c for 
(fci?a s ) _1 = —0.35, and the spin susceptibility is measured from the in situ imaging of 
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FIG. 7: (color online) Calculated T c , as a function of the interaction strength and the polarization 
rate P = [JV-j- — N±\/[Nf + N±\. In this figure, we assume the second-order phase transition 



dispersive speckle patterns. In the normal state above T c (the left side of the vertical line 
in Figj5]), Xetma agrees well with the observed spin susceptibility, without introducing any 
fitting parameter. 



While a good agreement with Ref. 47j is obtained, our result is somehow different 
from the experimental result done by Sommer and co-workers |53J]. In this experiment, 
the observed spin susceptibility in the normal state monotonically increases with decreasing 
the temperature. In contrast, the calculated spin susceptibility exhibits a peak structure, 
as shown in Fig. [6j This non-monotonic behavior is similar to the so-called spin gap 
phenomenon observed in the underdoped regime of high-T c cuprates 54J. In the present 
case, the decrease of x near T c is due to the development of the pseudogap in the single- 
particle density of states. For this discrepancy between the theory and experiment 53], 



55 



56| have recently 



although further analyses would be necessary, we note that Refs. 
pointed out that the experimental result may be understood by taking into account the 
non-equilibrium state associated with a quasi-repulsive interaction. 



IV. POLARIZED FERMI GAS IN THE BCS-BEC CROSSOVER REGIME 



We now consider the case of finite population imbalance. Figure [7] shows T c in the BCS- 
BEC crossover regime of a polarized Fermi gas (iV-f > N±), calculated within the framework 
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of ETMA. We briefly note that, since we are using the Thouless criterion in Eq. (fT2|) . 
the second-order phase transition is implicitly assumed. That is, possibility of the phase 
separation, which is accompanied by the first-order phase transition, is ignored in this figure. 

In the strong-coupling BEC limit, the system is well described by a mixture of N± tightly 
bound molecular bosons and — iVj_ excess "f^-spin atoms. Thus, the superfluid phase 
transition is dominated by the BEC of the former component. Since the phase transition 
temperature of an ideal Bose gas is proportional to N^ 3 (where is the number of bosons), 
T c in the extreme BEC limit is given by 



T c = T C BEC x 



2/3 

= T C BEC (1 - P) 2 / 3 , (17) 



.(N/2). 

where T C BEC = 0.218£f is T c in the BEC limit of a unpolarized Fermi gas. P = (iV-f — 
N±)/(Nf + iVj.) is the polarization rate. Equation ( fl7|) indicates that T c decreases with 
increasing P to vanish in the fully polarized limit (P — > 1). 

In the crossover region, as well as the BCS regime, FigJT] shows that T c vanishes at a 
certain value of P (= P c < 1). Since a polarized Fermi gas in the BCS regime is similar to 
metallic superconductivity under an external magnetic field, the vanishing T c at P c (< 1) is 
essentially the same as the suppression of the superconducting state by an external magnetic 
field. In the unitarity limit, one finds P c = 0.13, which is relatively close to the observed 



polarization rate P tc = 0.2 at the tricritical point of a 6 Li Fermi gas [58|, 159 ]. 

In the mean- field theory, the region of the phase separation (PS), which is surrounded 
by the first-order phase transition line, is obtained in the T — P phase diagram, as shown 
in Fig]8j (We summarize how to obtain this figure in Appendix B.) Since the mean-field 
theory is valid for the weak-coupling regime, the PS region would also appear in FigJTJ if 
one included the possibility of the first-order phase transition beyond the present treatment. 
To confirm this, however, we need to evaluate the thermodynamic potential Q, taking into 
account strong-coupling corrections within the framework of ETMA, which remains as our 
future problem. 

Figure ED^a) shows effects of the "effective magnetic field" h on the superfluid phase 
transition. As expected from the magnetic field effect on superconductivity, T c decreases 
with increasing h to vanish at a critical magnetic field h c . When we evaluate the polarization 
rate P along this T c -line, we obtain Figj9]^b). In this panel, T c is almost constant around 
h = 0, (dP / dh) h _> is close to the spin susceptibility x (> 0). Thus, in a sense, the positive 
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FIG. 8: (color online) Mean-field phase diagram of a polarized Fermi gas when (fcfa^ 1 = — 1. 
The solid line shows the second-order phase transition between the superfluid phase (SF) and 
the normal state (N). The dashed line shows the first-order phase transition, below which the 
phase separation (PS) of the superfluid phase and the normal state occurs. The tricritical point 
is obtained at Ptc = 0.199. When we ignore the PS phase and simply assume the second-order 
superfluid phase transition, we obtain the dotted line. In this case, one finds P c = 0.202. Assuming 
the second-order phase transition (solid line and dotted line), the reentrant region is obtain when 
0.164 < P < 0.202 (= P c ). 



P in panel (b) is a result of the correct treatment of the spin susceptibility in ETMA. 

In the inset of Figl^b), one sees a peak structure near the critical magnetic field h c . Since 
P c is given by this peak value, P c is found to obtain, not at h c , but below h c . As expected 
from the mean-field phase diagram shown in FigJHl one needs a more sophisticated treatment 
near P c and h c to include the first-order phase transition, as well as phase separation. 
However, apart from this, the origin of the peak seen in the inset of Figj9](b) is explained as 
follows. When the temperature T is fixed at a certain value, P monotonically increases with 
increasing h. On the other hand, when one decreases the temperature under the condition 
of a fixed h, the polarization P may decrease near T c , because of the suppression of the spin 
susceptibility due to the development of the pseudogap. (See Figj6j) In the case of FigJHl 
because the both mechanisms affect P, the polarization rate may decrease, when the latter 
effect becomes dominant. In particular, since the decreases of T c is most remarkable near 
h c (See Figfj^a).), this remarkable decrease of the temperature leads to the decrease of x, 
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FIG. 9: (color online) (a) Calculated T c , as a function of the "effective magnetic field" h. We take 
(fcfa s ) _1 = 0. (b) Polarization rate P at T c , as a function of h. The inset shows the polarization 
near the critical magnetic field h c at which the second-order phase transition disappears. The 
critical polarization rate P c is determined as the peak value seen in the inset. 

as well as P, as shown in the inset of Figj9]^b). We briefly note that, since the temperature 
is not fixed in panel (b), the negative value of (dP/dh)h~h c does not mean the violation of 
the required positivity of the spin susceptibility. In ETMA, the spin susceptibility is always 
positive, when the temperature is fixed. 

V. SUMMARY 

To summarize, we have investigated magnetic properties of a cold Fermi gas in the BCS- 
BEC crossover region. In the NSR theory, it is known that the spin susceptibility becomes 
negative in the crossover region. We showed that this unphysical result is also obtained in 
the ordinary (non-self-consistent) T-matrix approximation. We clarified that this negative 
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spin susceptibility originates from how to treat the self-energy correction and vertex correc- 
tion to the spin susceptibility. Improving this, we have succeeded in obtaining the positive 
spin susceptibility over the entire BCS-BEC crossover region. The calculated spin suscepti- 
bility agrees well with the recent experiment done by Sanner and co-workers 47], without 
introducing any fitting parameter. We have also applied our extended T-matrix theory to a 
polarized Fermi gas, and have calculated T c in the BCS-BEC crossover region. 

In this paper, we have considered the normal state above T c . Since the spin susceptibility 
is known to be strongly suppressed in the superfluid phase below T c , it is an interesting 
next challenge to extend the present theory to include the superfluid order parameter. This 
extension is also necessary in considering the first-order phase transition, as well as the phase 
separation, expected in polarized Fermi superfluids. 

In addition, we have only treated a uniform gas, for simplicity. Since a real cold Fermi 
gas is always trapped in a harmonic potential, the inclusion of this spatial inhomogeneity 
is crucial for detailed comparison of theoretical results with experimental data. However, 
since the key issues to overcome the negative susceptibility problem clarified in this paper 
is also valid for a trapped gas, our results would be useful for the further development of 
research for magnetic properties of trapped polarized Fermi gases. 
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Appendix A: Derivation of Eq. (1131) 

In the NSR theory, strong-coupling correction (= Xnsr) ^° the s P^ n susceptibility with 
0(S°) is given by 



flue _ d_ 
Xnsr - dh 



p,iu} n ,cr 



(Al) 

h-tO 
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Carrying out the /i-derivative, we obtain the contributions in FigsJU^b2) and (b3). Their 
expressions are given by, respectively, 



, <9 2 



/i=0 



-T E 



9^2 



n(q, ii/ n 



/i=0 



(A2) 



xSsr = 2T 2 E E r ^>^) K it (^„)] 2 K_ p>; (^ n -io; n )] : 



p,iw„ q.lVn 



h=0 



d 2 



IT E r (<?> Wn)x — a n(qr, ii/ n ) 



h=0 



(A3) 



In the BEC limit, the particle-particle vertex function in Eq. reduces to 60] 

r , • v 8tt 1 



(A4) 



Here, //b = 2/i + eb may be regarded as the chemical potential of molecular bosons, where 
€b = l/ma 2 s is the binding energy of a two-body bound molecule. Using the fact that the 
binding energy eb is very large in the BEC limit (a" 1 — » oo), one may expand <9 2 IT(qr, iv n )/dfi 2 
in Eq. ( 1A2I) with respect to e^ 1 . We then have 



dfi 2 



U(q,iu n ) = E 



d 2 f(£ P+q /2,t) 



W n - £p+q/2,t - ip-q/2,\. dlli 



+ 0(eZ 2 ) 



1 <9 2 A^ t 



+ 0(e b " 2 ), 



(A5) 



where = f(£, p ^) is the number of t-spin atoms in a free Fermi gas. Substituting Eq. 
fTA5l) into Eq. fTA2l . one obtains 



(b2) _ 2Td 2 N° 



Xnsr 



E T(q,iu n ). 



(A6) 



h=0 <J:^ n 



In particular, at T c , we find 



(b2 ) 16na s ( 2mT c 
Xnsr ~ m V 2tt 



'3^ 9 2 iV? 



2/ <9/i 2 



(A7) 



h=0 



We briefly note that, because d 2 U(q,ii> n ) /dfi^dfii is the order of e b 2 , one finds Xnsr 
0{e h 2 ). Thus, one can ignore Xnsr m ^ ne BEC regime. 
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Appendix B: Mean-field phase diagram of a polarized Fermi gas 



In the mean-field theory, the second-order phase transition is determined by solving the 
ordinary BCS gap equation at T c , 

u ^ WfaHM ; (Bi) 

together with the number equation, 

N = Y,f(^)- (B2) 

To evaluate the first-order phase transition temperature, we need to consider the ther- 
modynamic potential Q in the presence of phase separation (PS), which is given by 

fi(/i t ,/ij,,T, A,x) = xQsF(Pt,V±,T,A) + (1 - x)n N (^,ii h T). (B3) 

Here, Qsf and f^N are the thermodynamic potential in the superfluid (SF) phase and the 
normal state (N) region, respectively. Their mean-field expressions are given by 

fi ^-^ + Ek.-^41 -r5>(l + «-'W), (B4) 



p,cr 



fi N = -T^log(l + e-^ /T ), (B5) 



where E p>(7 = + A 2 — er/i is the Bogoliubov excitation energy. Since any intensive 
variable should have the same value in both the SF region and the N region in the PS phase, 
each of the chemical potential \x a and the temperature T takes the same value in f2sF and 
f^N- The superfluid order parameter A and the volume fraction x of the SF region are, 
respectively, determined from the stationary conditions of Q, 

° = 5A = ^f- < B6) 

= — = fi S F - On. (B7) 



Equation (IB 6ft gives the ordinary mean-field BCS gap equation. Equation (IB 7ft simply means 
Qsf = ^n- We solve Eqs. ( 1B6I) and ( 1B7I) . together with the number equations, 

JV t = xA^sf^t , // ; , A) + (1 - a;)iV tiN (/i t , /ij, (B8) 
N± = xN^spfa, m;, A) + (1 - s)^^, (B9) 
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to self-consistently determine A, x, \i a , below T c . In Eqs. (1B8j) and (IB 91) . A^sf and A^n are 
the number of tx-spin atoms in the superfluid region and the normal state region, respectively. 

The phase transition temperature T c from the PS phase to the normal state is obtained 
as the temperature at which the superfluid volume fraction x vanishes (x — 0). The phase 
boundary between the PS phase and the superfluid phase is determined by the condition 
x = 1. We have numerically evaluated these conditions to obtain the phase diagram in Fig. 
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